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ABSTRACT 



Aims. The EUV Imaging Spectrometer (EIS) instrument on board the Hinode satellite has access to some of the best coronal density 
diagnostics, and the high sensitivity of the instrument now allows electron number density, A'e, measurements to an unprecedented 
precision of up to ±5 % in active regions. This paper gives a thorough overview of data analysis issues for the best diagnostics of 
Fexii and Fexiii and assesses the accuracy of the measurements. 

Methods. Two density diagnostics each from Fexii (/1186.88//1195.12 and /1196.64//1195.12) and Fexra (.1196.54M202.04 and 
/1203.82//1202.04) are analysed in two active region datasets from 2007 May 3 and 6 that yield densities in the range 8.5 < 
log (A'e/cm"^) < 11.0. The densities are derived using v5.2 of the CHIANTI atomic database. Blending, line fitting, and instru- 
mental issues are discussed, and line fit parameters presented. 

Results. The Fe xn and Fe xm diagnostics show broadly the same trend in density across the active region, consistent with their simi- 
lar temperatures of formation. However, the high precision of the EIS measurements demonstrates significant discrepancies of up to 
0.5 dex in derived log A'e values, with Fe xn always giving higher densities than Fe xm. The discrepancies may partly be due to real 
physical differences between the emitting regions of the two plasmas, but the dominant factor lies in the atomic models of the two 
ions. Two specific problems are identified for Fexn /1196.64 and Fexin /i203.82: the former is found to be underestimated in strength 
by the CHIANTI atomic model, while the high-density limit of the /i203.82//i202.04 ratio appears to be inaccurate in the CHIANTI 
atomic model. The small grating tilt of the EIS instrument is found to be very significant when deriving densities from emission lines 
separated by more than a few angstroms. Revised wavelengths of 196.518 ± 0.003 A and 196.647 ± 0.003 A are suggested for the 
Fexin /1 196.54 and Fe xn /1 196.64 lines, respectively. 

Key words. Sun: corona - Sun: UV radiation 



1. Introduction 

Density is a fundamental plasma parameter that is accessible to 
astrophysical spectroscopists through the density diagnostic - a 
ratio of two emission lines that is sensitive to the electron density 
of the emitting plasma. High accuracy density measurements are 
of direct importance for constraining energy and pressure bal- 
ance within coronal structures, and also of indirect importance 
to the measure of the 'filling factor' of the plasma, i.e., the frac- 
tion (by volume) of the observed plasma structure that is actually 
emitting radiation. The latter is of significance when the spatial 
resolution of the telescope is not sufficient to resolve the physi- 
cal structures of interest, as is often the case even for our nearest 
star. The power of emission line density diagnostics lies in the 
fact that, when the plasma is optically thin, the derived density 
depends purely on the atomic data for the emitting ion and is 
free of any physical assumptions about the volume, abundances 
or temperature structure within the plasma. 

Density diagnostics for use in both analyses of stellar coro- 
nae and nebulae have been described and exploited for many 
ye ars. For the solar coro na, r eviews of early work are provided 
bv lDere & MasonI (1 19811) and lMariskal(ll992h . Steady improve- 



ments in the accuracy of atomic data calculations and access to 
high resolution, calibrated ultraviolet spectra from space-borne 
instruments have led to a number of detailed studies of density 
diagnostics from both individual ions and sequences of io ns, e.g., 
Brickhouse et al. ( 1995), Laming et al. ( 1997) and Youn g et all 
(119981) . 

A particular advance has been the routine access to most of 
the wavelength range 150-1600 A by the CDS and SUMER 
instruments on board the SOHO spacecraft for over a 12 
year period. A wide range of density diagnostics from many 
different ions has been exploited (e!g., feoschek e t al.L 1997 



Laming et al.L 119971; lYoung et al.L ll999HDel Zanna & Mason , 
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20031), however the density measurements from these two in- 
struments are generally accurate only to the 30-40 % level 
at best and often it is necessary to bin multiple spatial pix- 
els or exposures to obtain a useful signal. This is because the 
most commonly-used diagnostics generally involve only weak- 
to-medium strength emission lines and/or the diagnostics them- 
selves do not have high sensitivity to density. 

The best density diagnostics at typical coronal temperatures 
(1-3 million K) belong to the set of iron ions Feix-xv. They 
arise because of the complex atomic structures of these ions 
which lead to a number of metastable levels that are the source 
of the d ensity sensitivity of th e emission line ratios (e.g.. Sect. 
6.3.1 of lDere & Masonll 19811) . The density diagnostics of these 
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iron ions are mostly found below 300 A, a wavelength region not 
well covered by the SOHO spectrometers. 

The EUV Imaging Spectromer (EIS) on board the Hinode 
satellite - launched during 2006 September - observes two 
wavelength ranges below 300 A and, for the first time, scientists 
have routine access to high spectral and spatial resolution data in 
this region. A key factor is the high sensitivity of the instrument 
which is enabled b y a simple optical de sign that incorporates 
multilayer coatings jCulhane et al.Ll2007h . The sensitivity is par- 
ticularly high in the range 185-205 A where are found excellent 
density diagnostics from Fe xii and Fe xm. As demonstrated in 
the following sections, diagnostics from these ions yield density 
measurements with an unprecedented precision of up to 5 % in 
individual spatial pixels at exposure times of 30 s or less - a vast 
improvement over earlier instruments. 

The diagnostic potential of Fexii and Fexiii has been 
recognised for many years. The first so phisticat e d atom ic cal- 
culation for Fe XII was performed by iFlowerl (1 19771) and a 
number of UV dens i ty dia gnostic s hav e been studied (e.g., 
iFeldman et all Il983t iTavalet all 11991 . but for many years 
there were large discrepancies between predictions from atomic 
data mode ls and measured lin e intensities which were sum- 
marised in'Binello et al.' (2001'). Finally the large /^-matrix cal- 
culation performed by Storev et al. (2005) resolved many of the 
problems jDel Zanna & Mason! 120051) . The three emission lines 
used in th e present work (/1186 . 88, /11 95.12 and A196M) were 
shown by iDel Zanna & Maso ^ (120051) to give good agreement 
with other line ratios from Fe xii, but each line is affected by line 
blending. This will be discussed in detail during the course of 
the present woA^ 

For Fe xm, iFlower & Nussbaumej (1 197 41) calculated atomic 
data and produced line ratio calculations that yielded reason- 
ably consistent density measurements from a quiet Sun spec- 
trum. Detailed comparisons of the more recent atomic data 
models of Fe xiii in the CHIANTI atomic database (iDere et al.L 
Il997h with spe ctra from the SERTS rocket flight experiments 
( Neupert et al.L [T9 92) have been performed by lYoung et al.l 
0998^ and Landi (2002). In particular, the latter work found 
fliat the ^196.54//1202.04 and /1203.82//i202.04 ratios used in the 
present work give density estimates in good agreement with each 
other and also with other ions formed at similar temperatures. 

Early density results from EIS obtained using the Fe xii an d 
Fe XIII diagnostics hav e been presented in iBrooks et^ (l2007l). 
IWatanabe et all (l2007l). lOoschek et al.1 (l2007ah . I Doschek et al.l 
(l2007bl) and iDere et al.l (l2007h. In particular we highlight the 
density maps shown in Doschek et al.l ( l2007bl) that demonstrate 
the high quality of the EIS density measurements , and the den- 
sity comparisons shown in lWatanabe et al.l (|2007|) that show ex- 
cellent agreement between four different iron ion density diag- 
nostics. 

In this paper density measurements from two diagnostics 
each of the Fe xii and Fe xm ions are presented from two active 
region data sets. A careful discussion of line fitting and blending 
issues is given, and an estimate of the accuracy of the density 
measurements is made. A following paper will use the derived 
densities to estimate column depths and filling factors for the 
two ions. 



2. Instrument details 

The EUV Imaging Spectrometer (EIS) on Hinode observes 
two wavelength ranges: 170-212 A and 246-292 A. Th e in- 
strument is described in detail by iKorendvke et al.l (l2006l) and 



ICuIhane et al.l (l2007h . The key optical elements are: a parabolic 
primary mirror for focussing incoming radiation; a uniform line- 
spaced grating for dispersing the radiation and two back-thinned, 
EUV sensitive CCDs for detecting the radiation. Two aluminium 
filters lie on the optical path for blocking visible radiation, and 
a slit-slot assembly is situated between the mirror and grating 
containing four different-sized slits for selecting the spatial re- 
gion to be sent to the spectromer. The optical surfaces of both 
the mirror and grating are divided in two and coated with mul- 
tilayer coatings optimised for the two wavelength bands. Of the 
four entrance slits, two are available for detailed spectroscopy 
work with projected widths at the Sun of 1" and 2". Pixel sizes 
on the detector are 1" (spatial dimension) and 0.0223 A (spectral 
dimension). For single exposures the maximum field-of-view in 
Solar-Y is 512". 

For density measurements, the absolute radiometric calibra- 
tion of the EIS instrument is not important, only the relative cal- 
ibration between the wavelengths of the lines. The key uncer- 
tainty is whether the relative reflectivities from the mirror and 
grating between two wavelengths deviate from the values con- 
tained in the calibration pipeline. There are two components: 
the uncertainty due to variations in reflectivity across the optical 
surfaces; and a systematic error in the reflectivity curves due to 
some change in the instrument between the laboratory calibra- 
tion measurements and launch. For the former, measurements 
performed on flight-quality mirrors and gratings by See lv et aTl 
(12004 revealed that, when a beamline is scanned across the 
optical surfac es, the reflec tivity as a function of wavelength is 
found to vary. ISeelv et akl (j2004) averaged the resulting reflec- 
tivity curves and the standard deviations of the data points for 
the wavelengths considered here range from 3 to 8 %. In orbit 
the optical surfaces are fully illuminated which will serve to av- 
erage out the variations and so these standard deviation values 
will be over-estimates of the actual uncertainty in the relative 
calibration. In addition, these percentage errors are for the abso- 
lute reflectivity - the standard deviation for the relative reflec- 
tivities between two lines close in wavelength (e.g., the Fe xm 
/1202.04 and /1203.82) will likely be significantly less. For this 
reason we do not include any estimate of the relative calibration 
uncertainty in the error budget for derived densities. The error 
bars for the density measurements presented here are thus un- 
derestimates, but we believe the uncertainties in the relative cali- 
bration are comparable to or smaller than these errors and would 
not significantly impact the conclusions of this work. Of greater 
concern for the present work is whether there is a systematic er- 
ror in the shape of the multilayer reflectivity curves used in the 
EIS calibration. This can be investigated through the study of 
emission line ratios that are insensitive to the plasma conditions 
- an example being the Fe xii ratios considered in Appendix iBl 
The question of systematic errors is considered in the later sec- 
tions when the density results are presented. 

The basic unit of observation for EIS is a raster which is 
a consecutive set of exposures that either scan a spatial region 
(from solar-west to solar-east) on the Sun, or are fixed relative 
to the Sun's surface ("sit-and-stare" observations). A study is 
defined as a collection of one or more rasters that is designed 
to perform a particular science task. Each study and raster has a 
unique 20 character acronym that is used to identify it for science 
planning purposes. 

On account of data download volume restrictions, complete 
CCD images from each EIS exposure are not routinely transmit- 
ted to Earth. A key part of the raster definition process is thus to 
define 'wavelength windows' on the detector for selecting par- 
ticular emission lines and/or a reduced spatial coverage in the 
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Solar-Y direction. The widths of the windows need to be chosen 
carefully by the scientist so that accurate fitting of the emission 
line profiles will be possible. Typically window widths are cho- 
sen to be 24-32 pixels (0.54-0.71 A) but for some of the emis- 
sion lines considered here wider windows are necessary, and this 
will be discussed in the following sections. 



3. Atomic data 

To convert a measured line ratio value to a density value it is 
necessary to have a model for the atomic processes within the 
ion that predicts the relative strengths of the emission lines as a 
function of densit y. Here we use version 5.2 of t he CH IANTI 
atomic database dLandi et all 120061: iDere et al.L 1 19971) which 
contains atomic data and analysis software for calculating line 
ratios as a function of density. For Fe xii t he data co nsist of en- 
ergy le vels and radiative decay rates from 
( l2005h . electron collision strengths from 



Del Zanna & Masor 
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and proton collision rate coef ficients f rom Landman' ("1 9 781) . For 
Fe xm energy levels are from iPenn & Kuhn ( 1994), Jupe n et aTl 
([T993h and versio n 1.0 o f the NIST database; radiative de- 
cay r ates are from ^und (l2004i) . electron collision strengths 
from iGupta &Tavall (119981) and proton collision rate coeffi- 
cients fromlLandma^Q^ll). During the course of the present 
work, the CHIANTI team found an error in the CHIANTI 
electron collision stre ngth file for Fexm , whic h actually con- 
tained older data fromlFawcett & MasonI (Il989i) rather than the 
data o f iGupta & Tavall (Il998l) thaT were supposed to have been 
added in version 3 of the database jDere et al.L 12001 *). With the 
[Gupta & Taval (1998) dataset Fexiii yields significantly differ- 
ent densities (iLandill2002l) . 

The atomic transitions of the Fexii and Fexm density di- 
agnostic lines considered here are given in Table [1] and plots 
showing how the ratio s are predic t ed to v ary with density 
are shown in Fig. 3 of lYoung et alj (l2007bh. The diagnostics 
themselves are the three recommended by Young et"aLl (l2007bh 
- Fexii /1196.64/^195.12, and Fexm ^196.54/^202.04 and 
A203.82M202.04 - together with Fexii ^186.88^195. 12 which 
is comparab l e to /iI9 6.64//1195. 12but concern was expressed by 
lYoungetal.l(l2007bl) over the blending line S xi /1186.84. This is- 
sue will be discussed further in Sect. 15.11 Between them, these 
four density diagnostics are the best available in the EIS wave- 
length bands. The Fe xii diagnostics benefit from the high in- 
strument sensitivity at 195.12 A which leads to accurate mea- 
surements of the line ratios. The Fe xii diagnostics also benefit 
from a wide sensitivity to density from 10'^ to 10'^ cm"-' that 
arises through the shift of population from the ground '*5 term to 
first the term and then the term as density increases over 
this range. The EIS sensitivity for Fexm A2Q2.QA and /1203.82 
is lower than for the Fe xn lines (Table [TJ, but the value of the 
diagnostics lies in the strong sensitivity of the lines to density. 
E.g., the /1203.82M202.04 ratio varies by a factor 43 from 10'^ 
to 10'*' cm"-', whereas Fe xii /1 196.64//1 195. 12 varies by a factor 
8. This high sensitivity means that the error bars on a measured 
ratio translate to relatively small errors on the derived density. 

A more recent set of collisio n data for Fexm than th at in 
CHIANTI was published by lAggarwal & KeenanI (l2005b . and 
these have been compared with the v5.2 C HIANTI Fe xm mode l 
and SERTS rocket flight ob servations in Keenan et al. I (I2007h . 
The Fe xm model used by Keenan et al. (2007") consisted of 
collision data from Ag garwal & Keenan (2005) and radiative 
data from lAggarwal & KeenanI (l2004i) . The agreement with the 
CHIANTI model was found to be very good (although this 



was the CHIANTI model containing the Fawcett & Mason, 
1989, data rather than the model used here). Fig. [T] compares 
the theoretical ratio s from our revised CHIANTI model and 
iKeenan et al.l (l2007h for the two Fexm ratios. In bot h cases, for 
a given ratio value, the lAggarwal & KeenanI (|2005|) data yield 
lowe r densities by around 0.1-0.2 dex. The high-density limits 
from lKeenan et al.l (12007*) are higher than those from CHIANTI 
which has a significant impact on derived results from high den- 
sity pl asmas. These i s sues w ill be discussed further in Sects. |6] 
and |7] iKeenan et al.l (l2007l) prese nted high reso l ution spectra 
from the SERTS-95 rocket flight dBrosius et al.L 1 19981) where 
they were able to resolve the two components of the Fexm 
^203.82 blend, i.e., the lines at 203.79 and 203. 83 A. This is 
not possible with the lower resolution EIS spectra. IKeenan et"an 
(I2007D found that the /1203.79//1203.83 ratio was significantly 
diiJerent from theoretical predictions, with an observed value of 
0.20 and a theoretical value of 0.32. This raises concerns about 
the suitability of using the blended /1203.82 line as part of a den- 
sity diagnostic. We note, however, that summing the /1203.79 
and /1203.83 intensities from the SERTS-95 spectrum and form- 
ing a d ensity diagnos t ic wit h A2Q2.QA, the derived density (us- 
ing the IKeenan et al.L l2007l model) is log A^e = 9.1 which is 
consi stent with the other Fe xm densities presented in Table 9 
of Keenan et al.l ([20 07). We note the line widths measured by 
Kee nan et al.l (12007 ) for ^203.79 and ^203.83 are significantly 
different which is not expected from two lines of the same ion, 
and this may partly be responsible for the discrepancy between 
theory and observation for their intensity ratio. For the remainder 
of this work we will concentrate on the density results from the 
modified Fe xm CHIANTI model as CHIANTI is widely used 
amongst the solar physics community. 



-CHIANTI v5.2 
- Keenan et 
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Fig. 1. Theoretical variation of emission line ratios versus elec- 
tron density for Fexm /1203.82//1202.04 and /1196.54/^202.04. 
The solid lines show the results from the modified CHIANTI 
v5.2 mo del while the dashed lines show results from the Fe xm 
model of IKeenan et"aLl (l2007l) . 

The standard CHIANTI models for Fe xii and Fe xm include 
electron excitation/de-excitation, spontaneous radiative decay 
and proton excitation/de-excitation. Further processes that can 
be included are photo-excitation and stimulated emission which 
can be significant in low density conditions when electron col- 
lisions are less frequent. The method for including these lat - 
ter processes in CHIANTI is described in lYoung et alj (l2003h . 
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and examples of the e ffects of the processes on density d iag- 
nostics are provided in Flo wer & Pineau des ForetsI (Il973h and 
[Young & Esser ( 1999a b). For the solar case, the radiation field 
to be considered is the photospheric continuum which is mod- 
elled in CHIANTI as a black-body. For Fexii, such a radia- 
tion field has no effect on the density sensitive ratios consid- 
ered in this work. For Fexm there is no significant effect on 
A196M/A2Q2m, but there is an effect for A2Q3M/A2Q2M. 

Fig. |2] compares the predicted variation of Fexm 
/1203.82//1202.04 with density for two cases: when there is 
no incident radiation field (solid line) and when the Fe xiii emit- 
ting plasma is located 10" above the solar surface (dashed line). 
The radiation field serves to push the curve towards lower den- 
sities. A measured intensity ratio will thus yield a lower density 
when the radiation field is taken into account: for log Ng x 9 the 
reduction is about 0.05 dex. The change becomes less for higher 
densities, being negligible above log A^e = 10. The location 
of the Fe xiii emitting plasma above the photosphere is also a 
factor, with greater distances leading to a smaller effect. For the 
present work the effects of an incident radiation field will not be 
considered any further since the magnitude of the effect is small 
compared to some of the discrepancies discussed later, and also 
it is not possible to make a general statement about the location 
of the Fe xiii emitting plasma in the data sets considered. 

The present work is focussed on density diagnostics and 
so wavelengths and velocity measurements are not impor- 
tant. For all reference wavelengths quoted in the text, we 
use the the values from v.5.2 of the CHIANTI database 
dLandi et al.L l2006h. In particular we not e the revised Fexii 
wavelengths of iDel Zanna & MasonI ( l2005h which are included 
in CHIANTI and are valuable for interpreting the Fe xii /1195. 12 
line (Sect. 15. 2t . A by-product of the Gaussian fits to the emis- 
sion line profiles is line centroid measurements for each line, 
and these can be used to assess the accuracy of the CHIANTI 
wavelengths (Sect.jSJ. 
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Fig. 2. The theoretical variation of the Fexm /1203.82/.^202.04 
ratio as a function of density, as derived from CHIANTI. Two 
cases are shown: the curve calculated assuming no incident ra- 
diation field (solid line), and the curve calculated assuming the 
Fexm emitting plasma is located 10" above a photosphere of 
temperature 6000 K. 



4. Data calibration and error bars 

The EIS data presented here were all calibrated from the level-0 
FITS files using the EIS_PREP routine available in the Solarso^ 
IDL distribution. In addition to converting the measured CCD 
signal into calibrated intensity units, a key part of EIS_PREP is 
to flag bad data points. These can arise through pixel saturation, 
cosmic ray hits, or simply defective pixels on the CCD. In this 
work all such pixels were marked as 'missing' data and not in- 
cluded in the data analysis. The central outputs of EIS_PREP are 
two level- 1 FITS files, one containing calibrated intensities at 
each pixel, and the other containing error bars on these intensi- 
ties. The steps performed by the current version of EIS_PREP at 
the time of this analysis are described below in detail. 

The first step of EIS_PREP is to flag any saturated data. The 
EIS CCDs have a 14 bit dynamic range and so saturation oc- 
curs at 16,383 data numbers (DN). All such pixels are flagged as 
missing as described above. 

In the raw data, the spectra are found to sit on a background 
of around 500 DN that arises principally from the CCD bias, 
and secondly from the CCD dark current. It is not possible to 
estimate the CCD bias level directly for EIS data, so the bias 
and dark current levels are estimated directly from the science 
data as follows. For each 3D data window 2 % of the detector 
pixels are isolated that have the lowest DN values. The median 
DN value of these 2 % pixels is then set to be the background 
level and it is subtracted from the DN values of each pixel. 

Anomalously bright pixels are found on the EIS CCD images 
that arise from 'hot pixels', 'warm pixels' and cosmic rays. Both 
hot pixels and warm pixels are single pixels that have anoma- 
lously high DN values. A hot pixel is defined to be one that 
yields 25,000 electrons pixel"' s"' at room temperature (a spec- 
ification from the CCD manufacturer). Pixels that fall below this 
threshold but are still clearly identified as being anomalous when 
inspecting the data are referred to as 'warm' pixels. Maps of the 
locations of hot pixels are generated by the EIS team every 2- 
4 weeks following inspection of 100 s dark exposures and they 
are stored in Solarsoft. The hot pixel map that is closest in time 
to the science observation is used by EIS_PREP to mark the hot 
pixels as missing data. 

Before removing warm pixels and cosmic rays, the next step 
for EIS_PREP is to flag the pixels affected by dust on the CCD. 
Several small pieces of dust accumulated on the CCD before 
launch and are found to completely block the solar signal on the 
CCD at their locations. They are fixed in position and cover less 
than 0.1 % of the CCD, however two of the pieces do affect the 
strong lines Fe XI ^^188.23, 188.30 and Fe xii ^193.51 such that 
the lines can not be used over 15-30" spatial ranges in solar- Y. 

At the time of performing the present analysis, EIS_PREP 
did not specifically remove warm pixels, but many of these 
were removed naturally by the cosmic ray removal routine, 
EIS_DESPIKE. This latter is a wrapper routine that calls 
NEW_SPIKE, a routine developed for removing cosmic rays 
from 50j/0/Coronal D i agnostic Spectrorneter (CDS) data sets 
([Thompson et al.L 119981: iPike & Harrisoiil [2000). For CDS data 
processing it was typical for not only the identified CCD pix- 
els to be flagged, but also the nearest-neighbour pixels on the 
CCD. This is because there is often residual signal from the 
cosmic ray next to the brightest pixels. EIS sees significantly 
less cosmic rays than CDS apart from during the «5 minute 

' Solarsoft is a set of integrated software libraries, databases, 
and system utilities that provide a common programming and 
data analysis environment for Solar Physics. It is available at 
[http://www.lmsal.com/solarsofti 
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passes through the South Atlantic Anomaly, and the most use- 
ful function of EIS_DESPIKE is actually to flag wai'm pixels. 
Since warm pixels are only single pixel events, then the nearest- 
neighbour option is switched off for EIS. It is to be noted that the 
NEW_SPIKE routine was designed to be cautious when remov- 
ing cosmic rays from line profiles thus many weak warm pixels 
found within spectral lines are not removed, artificially enhanc- 
ing the emission line intensities at these locations. For the 2007 
May data sets analysed in the present work only around 2 % of 
the CCD pixels are warm pixels and so this is not a significant 
problem. 

The final step of EIS JPREP is to convert DN values into in- 
tensities in units erg cm"^ s"' sr"' A"'. The errors on the inten- 
sities are computed assuming photon statistics together with an 
error estimate of the dark current of 2.5 DN. 

4.1. CCD spatial offsets 

A further instrumental effect accounted for in the present work 
is the spatial off'sets between different wavelengths. Spatial off"- 
sets between images formed from lines on the two diff'erent EIS 
CCDs have been discussed by Young et al. (2007a). The emis- 
sion lines discussed here are exclusively from the short wave- 
length band and so these latter spatial offsets do not apply. 
However, the spectra on the short wavelength (SW) CCD are 
slightly tilted relative to the CCD's axes such that images from 
short wavelength lines are marginally higher on the CCD than 
long wavelength lines. The effect is due to a slight misalignment 
of the EIS grating relative to the CCD and which we refer to here 
as the "grating tilt". Note that this is not to be confused with the 
"slit tilt" (e.g..'Marisk a et afl.l2007h which is a misalignment of 
the EIS slits relative to the CCD. This latter effect results in slit 
images being slighly tilted on the detector, but does not result in 
image misalignment between lines at different wavelengths. 

The grating tilt is discussed in more detail in Appendix ICl 
where the gradient of the dispersion axis on the CCD is found to 
be -0.0792 pixels/A for the EIS SW band. Despite this small 
value, the effect on line ratios for lines separated by several 
angstroms can be significant as demonstrated in AppendixIC] For 
all of the density diagnostic lines considered here except Fe xii 
/1195.I2 the 2D intensity arrays that arise after Gaussian fitting 
has been performed have been adjusted to account for the spatial 
offsets. The intensities have been modified according to Eq. lC.ll 
and the intensity errors have also been modified accordingly. 



5. Line fitting - method and blending issues 

In the following sections we describe the methods used to 
automatically fit the density diagnostic emission lines with 
Gaussians. Firstly, though, we must ask two questions: is it nec- 
essary to fit Gaussians to the data? and is the Gaussian an appro- 
priate function to use? For density diagnostic work only the line 
intensity is needed, so it is simply necessary to sum up the inten- 
sity from the wavelength bins containing the emission lines with- 
out the need to resort to a model for the line profile. However, as 
will be seen, some of the lines are blended with other lines and 
so it is essential to have a model for the line profile in order to re- 
solve the blending components. Fitting the lines with Gaussians 
yields the line centroid and line width in addition to the intensity, 
which thus allow the density to be related to velocity and non- 
thermal broadening parameters of the solar plasma, providing a 
much more powerful data set for scientific analysis. The widths 



and velocities are also valuable for assessing the accuracy of the 
line fitting as will be demonstrated in Sect. [8] 

As to the appropriateness of Gaussian fits, in normal con- 
ditions the EIS line profiles are found to be very well repre- 
sented by a Gaussian function. In particular, no evidence is found 
for asymmetries or enhanced wings in the line spread function. 
However, non-Gaussian profiles due to real solar events are seen 
ofte n with E IS and examples are given in Imada et al. (2007) and 
Chi for et al.l ([2008). Typically such events are isolated and can 
be readily identified following an automatic Gaussian fit to the 
data, as they will have anomalous line widths and/or velocities. 
For the data analysed in the present work there is no evidence 
for any significant deviations from Gaussian profiles. 

The high quality of the EIS data mean that the lines con- 
sidered here can be automatically fit with Gaussians at virtually 
every spatial pixel in the EIS images. However, care has to be 
taken in applying automatic fitting routines due to the influence 
of both nearby and blending emission lines. In this section we 
go through each line and highlight blending and fitting issues. 
A key point is that the fits must be accurate and robust for the 
wide range of conditions encountered in the solar atmosphere. 
The discussion in thi s section bui l ds on t he initial overview of 
the lines presented in lYoung et al.l (l2007bl) . 

To fit Gaussians to the EIS Unes the MPFIT procedures of 
C. MarkwardQ have been implemented to automatically fit the 
calibrated line intensities from the spectral line windows out- 
put by the IDE routine EIS.GETWINDATA. The spectrum from 
each image pixel is fitted, with data points weighted with the 
Icr error bars calculated by EIS_PREP. The routine developed 
by PR. Young for fitting single Gaussians to the spectral lines 
is called EIS_AUTO_FIT and is available in the Solarsoft IDE 
distribution. For some of the spectral features discussed below 
it was necessary to perform multiple Gaussian fits and custom 
versions of EIS_AUTO_FIT were created for these. 

In the sections that follow the reference wavelengths of emis- 
sion lines (denoted by, e.g., /1195.12) are taken from v5.2 of the 
CHIANTI atomic database. Wavelengths are given to two signif- 
icant figures, except where greater accuracy is required. Two of 
the density diagnostic lines are actually self-blends of two tran- 
sitions from the ions and we use /1186.88 to denote the blend of 
Fexii /1186.854 and ^186.887, and ^203.82 to denote the blend 
of Fe XIII Aim. 191 and ^203.828. 

5.1. FeXI1 186.88 

The emission line at 186.88 A is generally strong and well- 
resolved in active regions (Fig. |3]l and it largely comprises of 
two Fexii transitions with CHIANTI wavelengths of 186.854 
and 186.887 A. Sxi /1186.839 is a known blend and its contri- 
bution can be assessed from the /il88.617 or /1191.266 transi- 
tions which are related by branching ratios as they arise from 
the same upper level (Islp^ ^Si). From CHIANTI the ^186.84, 
/il88.68 and A191.21 intensities lie in the ratios 0.195:0.580:1.0. 
In the data sets discussed in this work /1 19 1.27 is not observed, 
but /1188.68 lies in the window containing the strong Fe xi lines 
at 188.22 and 188.30 A (Fig.©. 

To demonstrate the contribution of S xi to the Fe xii /1186.88 
fine the intensity of S xi /1188.68 in the May 3 and May 6 data 
sets considered in Sects. |6]and|7]has been measured in the Fe xi 
/1188.22 window. Fig.|4]shows the /1188.22 wavelength window, 
where a total of six emission lines can be seen. A six Gaussian fit 
was applied and, to ensure robust fits, each line was assumed to 

" http://astrog.physics.wisc.edu/~craigm/idl/idl.htmI. 
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have the same width. Otherwise ampUtudes and centroids were 
free to vary. In both the May 3 and May 6 active region datasets, 
the six Gaussian fits were generally accurate. A 35 pixel high re- 
gion on the detector is affected by the presence of a dust particle 
in the May 6 data set and this prevents accurate fits to the six 
lines in this region - hence the gap in the data points in the right 
panel of Fig.|5] There are also isolated pixels where poor fits are 
found, typically due to missing pixels or weak lines. 

With the intensity of S xi /il88.68 measured, one can esti- 
mate the intensity of S xi /il86.84 by simply multiplying by the 
factor 0.34. This derived intensity can then be compared with 
the measured Sxi-Fexii feature at 186.88 A. For the two data 
columns in the May 3 and May 6 datasets analysed in Sects. |6] 
and 111 the results are shown in Fig. |5] The key point to note is 
that the majority of the pixels show a S xi contribution of < 5 %. 
Where the electron density is high in the May 3 data (pixels 50- 
170, see Fig.fTvTi the S xi contribution is around 2 % as the Fe xii 
lines are more sensitive to higher densities than the S xi line. 

The contribution of the S xi line at the 5 % level would have 
the most effect on the Fe xii densities where the Fe xii ratio is 
most sensitive, i.e., around log A^e = 10.5. Using CHIANTI we 
find the Fe xii density is reduced by 0.08 dex. At l og A^e = 9.0 
the den sity is reduced by 0.05 dex. We also note that Brown et al.l 
(120081) list a Nixiv Hne as possibly blending with S xi Al 88.68, 
potentially compromising the use of the S xi line for estimat- 
ing the /1186.84 contribution. Together with the fact that Sxi 
/1188.68 is a weak line that is difficult to measure, we have de- 
cided not to correct Fe xii /1186.88 for the S xi contribution in the 
present analysis. 
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Fig. 3. An example double-Gaussian fit to the Fe viii /II 86.60 and 
Fe XII /1186.88 lines from the May 6 data set. The chosen spatial 
pixel is indicated. The data are plotted as a thin line and the 
fit with a thick line. The wavelength pixels used in the fitting 
are indicated by crosses. Note that the weak line on the long 
wavelength side of Fexii /1186.88 has not been included in the 
fit. 



For fitting the 186.88 A feature, it is necessary to consider 
nearby lines that are not directly blended but can affect line fit- 
ting. Fe vm /1 186.60 is not a blend, but sufficiently close that it 
can affect background estimates when fitting /1186.88. For this 
work both lines have been fitted simultaneously with two inde- 
pendent Gaussians. A weak line is found in the long wavelength 
wing of /1186.88 at 186.98 A which has been suggested to be a 
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Fig. 4. Spectrum from a single spatial pixel in the May 3 
dataset showing the wavelength window used to obtain the Fe xi 
/1/1188.22, 188.30 feature. A six Gaussian fit to all of the lines in 
the window is shown. The weak S xi /1188.68 line used to esti- 
mate the contribution of S xi to Fe xii /1186.88 is indicated. 



Ni XI transition jBrown et al 1 I2008I) . Checks on the May 3 and 
May 6 data sets have demonstrated that this line is always weak 
relative to /1186.88, at most around 20 % of the intensity. Adding 
an extra Gaussian to the fit to account for this line often works 
well, however there are a number of pixels that yield poor fits 
if there are missing data or weak lines. For the present work it 
was decided to ignore the Ni xi line by omitting the pixels at this 
wavelength - an example can be seen in Fig. [3] This method 
leads to much more robust fits while still yielding accurate esti- 
mates of the /1186.88 intensity. 

5.2. FeXI1 195.12 

The Fe xn emission line at 195. 12 A is usually the strongest Une 
observed by BIS and it dominates the spectrum in the region 
194-196 A. Despite this, some care has to be taken when auto- 
matically fitting the line. Firstly, there is another Fe xn transition 
at 195.18 A which lies within the /1 195. 12 profile and thus dis- 
torts the observed line profile from a Gaussian shape. This line is 
predicted by CHIANTI to be < 10 % of the strength of A 195. 12 
below 10'" cm"^, but this increases to 22 % at 10" cm"^ and so 
can not be ignored if the data set shows high densities. A detailed 
discussion of the effects of the blending /1 195. 18 line is provided 
in Appendix lAl 

In addition to the Fe xii blending component there are a num- 
ber of weak lines either side of /1 195. 12 in the spectrum, as can 
be seen in the two panels in Fig. |6] Simply performing a single 
Gaussian fit over the wavelength range 194.80-195.50 A, for ex- 
ample, will lead to an incorrrect fit to the spectrum background 
which will result in the intensity of the /1 195. 12 Une being un- 
derestimated at the 1-3 % level. 

The solution to these two problems in the present work is 
as follows. Spectral regions are identified that are free from 
emission lines and are used to constrain the continuum level. 
The regions containing the weak lines are not included in the 
fit. This is illustrated in Fig. |6] where crosses indicate the pix- 
els that have been used for the Gaussian fitting. To deal with 
the /1 195. 18 blend a two Gaussian fit has been performed. The 
stronger /1 195. 12 line has free parameters for the centroid, width 
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and amplitude, while /! 195. 18 is forced to be 0.06 A to the long 
wavelength side of /II 95. 12 (the CHIANTI wavelengths for the 
two transitions are 195.119 and 195. 179 A), and to have the same 
width as /1 195. 12. The amplitude of /1 195. 18 is free to vary. 

While /1 195. 12 is the strongest Fexii line observed by EIS 
and thus likely to be the most commonly used line for density 
diagnostics, there are two other strong Fexii lines at 192.39 A 
and 193.51 A that are density insensitive relative to /1 1 95 . 1 2 and 
thus can be substituted for this line to be used as density diagnos- 
tics relative to /1186.88 and /1196.64. These lines are unblended 
and are less likely to saturate on the detector than /1 195. 12 since 
they are both weaker transitions and the EIS effective area is 
lower at these wavelengths. /1192.39 and /1193.51 are discussed 
in Appendix |B] where it is shown that the ratios of /1192.39, 
1 93 .5 1 and /1 1 95 . 1 2 closely match the predictions of CHIANTI 
and thus they all can be used in density diagnostics. 
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Fig. 6. Two example spectra from the May 3 data set showing 
the emission lines neighbouring the strong Fe xn /1195. 12 line. A 
two Gaussian fit has been performed to /i 195. 12 and the neigh- 
bouring Fexii /1195.18 line. The crosses show the wavelength 
pixels that have been included in the fit, and the thick line shows 
the fit result. A thin line is used to highlight the weaker Gaussian 
that represents /1195.18. Note that there are a number of weak 
lines neighbouring the Fe xii feature that potentially affect the 
fits. The peak of the /1195.12 line extends by a factor five be- 
yond the top of each of the two plots. 



/1196.64 there are three lines: S x /t l96.81, Fexii /II 96 92 and an 
unknown line at around 197.02 A dBrown et all l2008 l). Beyond 
this latter line there is a small region at around 197.2 A which 
is free of lines and suitable for estimating the background level. 
On the short wavelength side of /il96.54 there is a Fe vii line at 
196.42 A dBrown et al.L 2008) but this is generally negligible in 
active region conditions. 

For the fits performed in this work, only the pixels directly 
neighbouring the Fe xii and Fe xiii lines were used for estimating 
the background level (Fig. IT) as the wavelength window was not 
wide enough to observe the background region around 197.2 A. 
A uniform background was assumed for the fitting rather than a 
linear fit in order to prevent the S x /il96.81 line distorting the 
derived background. 

In regions where there is an absence of hot plasma and the 
density is low, the Fexiii /1196.54 line can be very weak and 
thus difficult to fit, but this was a not a problem for the data sets 
considered here. 

iBrown et all (|2008) list a Fe viii line at 196.65 A which thus 
blends with Fe xn /i 196.64. CHIANTI predicts the Fe viii line to 
be a factor 0.06 weaker than Fe viii /1194.66 (visible in Fig.|6]l, 
and thus it would be expected to be negligible in most cir- 
cumstances. However, in active regions Fe vm is often seen to 
be strongly enhance d in small brightenings or loop footpoints 
dYoung et al.L l200'7ah . and in these regions /I196.65 can be sig- 
nificant. Examples are found in both the May 3 and May 6 data 
sets where methods of correcting the measured /1196.64 intensity 
for the Fe viii contribution are discussed. 
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Fig. 7. An example double-Gaussian fit to Fexiii /1196.54 and 
Fe XII /1196.64 from the May 6 data set. The chosen spatial pixel 
is indicated. The data are plotted as a thin line and the fit with a 
thick line. The error bars on the intensity measurements are also 
shown. 



5.3. FeXIII 196.54, FeXII 196.64 

These two lines are fitted simultaneously with a double 
Gaussian, with each parameter free to vary. The main difficulty is 
to identify a clean continuum area in the spectrum for accurately 
fitting the background for the fits. To the long wavelength side of 



5.4. FeXIII 202.04 

Fexiii /1202.04 is unblended and can be fit with a single 
Gaussian. We choose to fit a single Gaussian to the wavelength 
region 201.90-202.32 A (Fig.®. 
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Fig. 8. An example single Gaussian fit to the Fe xm /1202.0 emis- 
sion Une. The May 3 data set is used, and the image pixel indi- 
cated. 



5.5. Fe XI 1 1 203.82 

This feature is a self-blend of two Fe xiii lines with CHIANTI 
wavelengths of 203.797 and 203.828 A. The two lines have weak 
density dependence (Fig. [TOb . with the ratio almost constant 
above 10^^ cm"^. There is an additional line due to Fe xii found 
at 203.728 A which is not clearly separated from the Fe xm lines 
(Fig.|9]l. Generally this line is always weaker than the Fe xm fea- 
ture, but in quiet Sun conditions it can be comparable in strength. 

The fitting method chosen for these lines is to fit the wave- 
length region 203.5-204.1 A with three Gaussians each set to 
the same width (which is free to vary). The two Fe xm compo- 
nents are assumed to have a fixed ratio /1203.797//1203.828 of 
0.40. Although Fig. [TOl demonstrates that this is only valid for 
log Ne; > 9.0, the assumption is preferable to allowing the line 
ratio to freely vary when performing automatic fits to the data. 
The close separation of the two Fe xm lines also means that their 
precise ratio is not crucial to fitting the spectral feature. The peak 
of the Fe xii component is free to vary in the fit, but the centroid 
is fixed relative to the 203.828 A line using the separation from 
the CHIANTI database. Fe xm /1203.797 also has a fixed centroid 
relative to the /i203.828 line. In all there are four free parameters, 
with two additional parameters for a linear fit to the background. 



6. Quiescent active region data set (2007 IVlay 6) 

The data sets presented in this and the next section are from 
active region NOAA AR 10963 which rotated across the so- 
lar disc during 2007 April 24 to May 9. On May 6 the re- 
gion was close to the west limb and observed by EIS with the 
study AR_VELOCITY_MAP which consisted of an initial con- 
text raster with the 40" slit, followed by a long duration 1 " slit 
raster covering a region 330" x 304". Fig. [TT] shows an image 
from the 1 " raster formed from Fe xn /1195. 12. The vertical black 
bands arise from Hinode passing behind the Earth as the raster 
progressed (the observation was taken during the Hinode eclipse 
season). Gaussian fits were performed on the spatial region high- 
lighted in Fig.[TT] The region has size 21" x 304" and it crosses 
a filament and some coronal loops. At each pixel in this sub- 
region each of the density diagnostic lines discussed in Sect. |5] 
has been fitted. We consider a slice through the data set at X- 



Fig. 9. An example fit to the feature at 203.8 A that comprises 
Fe XII /1203.72, and Fexm ^^203.79, 203.83 from the May 6 data 
set. The spatial pixel chosen is indicated. Three thin lines are 
used to show the three individual Gaussians from the fit. 
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Fig. 10. Theoretical variation of the Fe xm /1203. 80/^203.8 3 ra- 
tio derived from CHIANTI. 



pixel 10 and present the derived densities from the Fexii and 
Fe XIII fi nes in Figs.[T2land[T3l respectively. Densities calculated 
with the Keenan et al. 1 (120071)" Fe xm atomic model are shown in 
Fig. US] (available online). The upper panel in each plot shows 
the images derived from the Fe xii /1195.12 and Fexm /1202.04 
lines, respectively, with the selected image column indicated. 

Considering first the Fexii ratios /1186.88//1195.12 and 
/1196.64//il95.12, the line intensities from the fitting methods 
described in the previous sections are used for each line. /1186.88 
is not corrected for the Sxi /tl86.84 contribution (Sect. 15.1b . 
but /1196.64 is corrected for the Feviii /1196.65 contribution 
(Sect. 15.3b . The correction was performed by fitting the Fe viii 
/1194.66 emission line at each pixel in the data set. The fitting 
method for this line was similar to the nearby Fexn /1 195. 12 
(Sect. 15.2b . with only pixel sections chosen for the line and 
the background regions at 194.5, 195.7 and 195.9 A. From 
CHIANTI, the Feviii ^196.65//1194.66 ratio is insensitive to 
density and has a theoretical value of 0.058. Thus the Fexii 
/1196.64 intensity was derived by subtracting 0.058 times the 
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Fig. 11. Image of AR 10963 obtained in the Fe xii ^195. 12 line 
on May 6. The raster began at 18:19 UT and completed at 23: 13 
and was obtained during the Hinode eclipse season. The three 
thick, black, vertical lines show when the satellite passed behind 
the Earth. The two thin vertical black lines indicate the spatial 
region used for the analysis in Sect.|6] 

intensity of Fevm /1194.66 at each pixel. The maximum con- 
tribution of Fe VIII to the measured line intensity of the 196.64 A 
feature is found to be 9 % using this method. 

The measured intensity ratios were converted to densities 
using the CHIANTI database, and are plotted in Figs. [12] and 
[T3] The Icr error bars on the intensities have been translated 
to errors on the derived densities, and thus vertical lines are 
used on these plots, connecting the lower limit to the upper 
limit. The densities track each other reasonably well along the 
slit, with the positions of peaks and troughs in good agreement. 
However, the /1196.64//1195.12 ratio generally gives higher den- 
sities than /1186.88//1195.12, and particularly so at certain of the 
peaks along the slit, e.g., around pixels 170, 200 and 240. 

The dashed line in Fig. [12] shows the variation of Fe vm 
/1194.66 along the slit, and it is clear that the largest discrepan- 
cies between the two Fe xii ratios occur where Fe vm /1194.66 is 
strongest, implying that the Fe vm contribution to Fe xii /1196.64 
is not fully accounted for. This suggests that there are problems 
with the Fe viii atomic model in CHIANTI. It is to be noted that 
adjusting the theoretical Fe vm /1196.65//1194.66 ratio to around 
0.25 brings the ^196.64/^195. 12 and ^186.88//1195.12 densities 
into good agreement. 

The two Fe xiii density sensitive ratios show excellent agree- 
ment (Fig. [T3l l. except in the pixel region to 30 where the 
/1203. 82/^202. 04 ratio is lower by 0.2-0.3 dex, although the er- 
rors bars on the /1196.54//1202.04 ratio are significantly larger 
in this region. At all pixels the /1196.54//1202.04 ratio yields a 
marginally higher density than /1203.82//1202.04. Densities de- 
rived from the Keenanet al. (2007) atomic model (Fig. [TSll are 
very similar, being around 0. 1 dex lower than for the CHIANTI 
model, and the two ratios are in even better agreement, yielding 
almost identical densities along the EIS slit. 

In Fig. [14] we compare the densities from Fexii 
/1186.88M195.12 and Fexra /1203.82M202.04. (We choose 
these two ratios for the comparison due to the Fe vm correction 
problem for /1196.64//1195.12 discussed earlier, and the fact that 



Fexiii /1203.82 is a stronger line than /1196.54 at low densities 
and thus yields densities with smaller error bars.) The Fexii 
ratio yields higher densities at all pixels by up to 0.4 dex, 
and we also note that Fexiii generally shows less variation 
pixel-to-pixel than Fexii. These differences will be discussed 
further in Sect. [9] 

7. High-density data set (2007 IVIay 3) 
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Fig. 16. Context image of AR 10963 obtained on 2007 May 3 
00:03 UT, shortly before the raster analysed in Sect. [7] The im- 
age was obtained by making a six exposure raster with the 40" 
slit. The image is formed in the emission line Si x /1261.06 which 
is formed at a similar temperature to Fe xii /1 195. 12. The box in- 
dicates the location of the narrow slit raster that followed the 
context raster. 

The May 6 data set demonstrates the use of the Fe xii and 
Fe XIII diagnostics in a quiescent part of an active region, typ- 
ical of large, developed coronal loops. The two sets of ion ra- 
tios are sensitive to higher densities, and so the EIS data archive 
was searched for examples of high densities. On 2007 May 3, 
the same active region showed small flaring activity and this 
was captured with the EIS study DIAG_40X180_S0_30S which 
consists of two rasters, the first being the same 40" slot context 
raster used by the AR_VELOCITY_MAP study, and the second 
being a raster covering an area 40" x 176" using the 1" slit, the 
latter having an identical line list to the narrow slit raster used on 
May 6. DIAG_40X180.S0_30S was run during the period 00:07 
to 05:28 UT with the slot context raster run once at the start of 
the period, and the narrow slit raster repeated 15 times. A GOES 
class C9 flare began at May 2 23:15 UT, peaked at 23:50 UT 
and decayed by May 3 02:00 UT. The first narrow slit raster ran 
during the period 00: 1 1 UT to 00:32 UT during the decay phase 
of the flare and the flare line Caxvii /1192.82 is prominent in 
the spectrum. Fig.[T6]shows an image obtained with the context 
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slot raster at 00:07 UT and reveals the large-scale structure of 
the active region. The image is from Six /1261.06 as the Fexii 
and Fexiii lines are saturated due to the flare in this data set. 
Saturation of /1195.12 also occurs in the brightest parts of the 
narrow slit raster image, and so for our density study a data col- 
umn is selected that is to one side of the very brightest regions 
in the image (indicated in the upper panels of Figs. [TT] and [TSll. 
Even in this data column, one or two pixels in the centre of the 
/1195.12 profile are saturated at the brightest part of the image, 
over Y-pixels 54 to 69. Surprisingly, though, the two Gaussian fit 
to the 195 A feature still yielded good results. This was checked 
in two ways. Firstly the insensitive /1192.39//1195.12 ratio was 
plotted along the data column (Appendix iBl) and no significant 
feature is found over the saturated region. Secondly, the density 
was re-computed using the /1186.88//1192.39 ratio, and differ- 
ences of at most 0. 15 dex were found. 

As for the May 6 data set, Fe xii /II 96. 64 has been corrected 
for the Fe viii blend using the /1194.66 line, although the Fe viii 
contribution to the /1196.64 feature is at most 5 % here. 

X-pixel 10 in the May 3 data-set was selected for the den- 
sity analysis, and the densities from the pairs of Fe xii and Fe xm 
ratios are pres ented in Figs. [ T7] and [TS] The Fexm densities de- 
rived from the iKeenan et al.l (|2007|) atomic model are presented 
in Fi g . l20l (available online). 

The derived Fe xii densities again track each other reasonably 
well, with /1196.64//1195.12 yielding higher densities by up to 
0.4 dex, with the largest discrepancies being at high densities 
(Fig.fTTIl. A strong peak in the Fe viii /1194.66 line is seen around 
Y-pixel 57 which may partly account for the density discrepancy 
here, but at other locations this can not be the case. 

For Fe xm the two ratios again track each other very well at 
low densities (Fig.[T8ll, but at high densities /l203.83//i202.04 is 
found to get close to the high-density limit of log A^e - 11-3 
of the ratio (indicated by the large error bars). This does not 
happen for /1196.54//1202.04 which has greater sensitivity to 
high densities. This issue will be discuss ed further in Sect. [9l 
The same problem occurs when using the IKeenan et alj ( |2007|) 
atomic model (Fig.|20l), although as the .1203. 82/.1202. 04 high- 
density limit is slightly larger for this model (Fig. [1]) then the 
discrepancy between the two ratios is not as large at high den- 
sities. Densities from /1196.54//1202.04 are 0.2 dex lower at the 
highest densities using the Keenan et al. (2007) model compared 
to the CHIANTI model. 

In Fig. [19] the densities from Fexii .1186.88M195.12 
and Fexm ^196.54//1202.04 are compared. .11 96. 54/^202. 04 
is preferred to /1203.82//1202.04 for Fexm as we believe 
/1203.82//1202.04 is not accurately measuring the density for 
log A^e ^ 10.0 (Sect.|9]l. The discrepancies noted for the May 6 
data set are also seen h ere, with Fe xn yieldi ng higher densities 
by up to 0.5 dex. If the IKeenan et alj (l2007h model is used the 
difference between Fe xii and Fe xm rises to 0.7 dex. The magni- 
tude of this difference wiU be discussed in Sect.jH 



8. Integrity of the line fit parameters 

The density analysis presented in the previous sections used the 
emission line intensities that are derived from Gaussian fits to 
the emission lines. The fits also yield centroid and line width 
information for each line and these are valuable for assessing 
the quality of the line fits and for revealing any unaccounted for 
line blending. Each of the three Fe xii lines and each of the three 
Fe xm lines would be expected to show the same line width and 
velocity behaviour across the data and the plots in Fig.l^Tlshow 



comparisons for the various Une pairs. For Fexii, the compar- 
isons are done against /1 195. 12 (the eight panels on the left side 
of Fig. I2TI1. while for Fexm they are done against /1202.04 (the 
eight panels on the right side of Fig. |2TI) . It is to be noted that 
the measured /1186.88 line is a blend of two Fe xii emission lines 
separated by 0.033 A which thus affects the width and veloc- 
ity measurements. Fexm /1203.82 was assumed to be a blend 
of two Fe xm lines when calculating densities, however, in the 
Gaussian fitting (Sect. 15.51 ) the line properties of both compo- 
nents were derived thus the comparisons in Fig. |2T|are for the 
longer wavelength /1203.828 component. 

Considering first the line widths (the full width at half maxi- 
mum, FWHM, in this case) we plot the widths of the line of in- 
terest against the reference line in row 1 and row 3 of Fig.l^Tlfor 
the May 6 and May 3 data sets, respectively. Ideally the points 
should lie on a diagonal line running from the bottom-left to 
top-right, indicating that all the Fexii or Fexm lines broaden 
together For the May 6 data set the line width of /1 195. 12 is 
very consistent along the slit, with an average width of 66.0 + 
1.3 mA. However, both /1186.88 and /1196.64 show a signif- 
icantly broader spread, with average values of 78.0 ±4.1 mA 
and 67.3 + 4.5 mA. The greater spread is simply due to the lower 
signal-to-noise of the two lines. The /1196.64 and /1195.12 widths 
are in excellent agreement, while the larger width of /1186.88 
arises from the fact that it is a blend of two lines with wave- 
lengths 186.85 and 186.89 A. 

For Fexm, the greater spread of widths for /1203.83 and 
/1196.54 are again explained by the lower signal-to-noise in these 
lines. The average width for /1202.04 is 67.4 +1.7 mA, while for 
.1203.83 and .1196.54 they are 69.9 + 3.4 mA and 63.2 + 5.0 mA, 
respectively. The latter is noticeably lower than the other width 
measurements, although it is consistent within the error bars. 

Moving to the May 3 data set, a greater spread of line widths 
is found particularly for Fe xii. The average widths for .1195.12, 
.1186.88 and .1196.64 are 67.9 + 3.1, 84.0+4.4 and 69.4+4.2 mA, 
respectively. The average .1186.88 width is significantly larger 
than for May 6 due to an increasing contribution of the 186.85 A 
component at high densities. The other two lines have compara- 
ble widths to May 6, only slightly broader as might be expected 
for this flare data set. The average widths for Fexm .1202.04, 
.1196.54 and .1203.83 are 67.5 ±2.8, 66.0+3.9 and 71.0±3.1 mA, 
respectively. We note that, as with the May 6 data set, .1203.82 
is a little broader than the other lines which implies the three 
Gaussian model for the lines at this wavelength may not be cor- 
rect. Either an unaccounted for blending line, or errors in one or 
more of the rest wavelengths could be responsible. The effect, 
however, is small. 

For comparing the line centroid measurements we consider 
the line velocities as given by 



where c is the speed of light, A is the measured emission line 
centroid and .Iret is the reference wavelength of the line (taken 
from CHIANTI). The positions of the EIS emission lines on th e 
detector vary over the orbit by ^ ±1 pixels dBrown et al.Ll2007l) . 
and so the measured centroids are normalised here relative to the 
reference lines (.1195.12 for Fexii and .1202.04 for Fexm), i.e., 
the wavelengths of these lines along the EIS slit are assumed to 
average to the reference wavelengths of the lines. As for the line 
widths, the line velocities should ideally fall on a diagonal line 
running from bottom-left to top-right in the figures, indicating 
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that each hne should reveal the same blueshift or redshift. In ad- 
dition the diagonal line should pass through (0,0) if the reference 
wavelengths of the lines are both correct. 

For the May 6 data set there is little dynamic activity along 
the slit which is revealed in the /1202.04 and /1195.12 measure- 
ments which each have a very narrow spread of values. The other 
lines show a broader spread of values due to the lower signal- 
to-noise of these lines. Fexii /1196.64 shows a small clump of 
values around -h20 km s"' which is caused by the blending 
Fe vm line mentioned in Sect. 15.31 Of the four comparison plots 
for May 6, only for /1203.83 are the points close to the (0,0) 
point, indicating accurate reference wavelengths. For /1186.88 
the ^ -20 km s ' offset is explained by the blending 186.85 A 
line. For Al 96.64 the average offset is 10.9 + 3.4 km s"' which 
suggests a revised reference wavelength for the line of 196.647 A 
(assuming the CHIANTI wavelength for /1195.12 is correct). 
This compa res with the value o f 196.645 A in the quiet Sun 
spectrum of iBrown et al. I (I2008h . The average offset for Fe xm 
/1196.54 is -33.9 + 3.9 km s'', suggesting a revised reference 
wavelength of 196.518 A (assuming the CHIANTI wavelength 
for /1202.04 is correct). The average value in the five spectra pre- 
sented bv lBrown et al.l ( l2008l) is 196.520 A. The uncertainties in 
the reference wavelengths of Fexii /1195.12 and Fexm /1202.04 
are +0.002 A (iDel Zanna & Masonl l2005t iBrown et all l2007l) . 
thus our revised wavelengths for the Fexii /1196.64 and Fexm 
/1196.54 are 196.647 + 0.003 A and 196.518 + 0.003 A, respec- 
tively. 

The line velocity plots from the May 3 data set in the bottom 
row of Fig.|2T]are quite striking. There is clearly significant dy- 
namic activity in the data set which is actually related to the high 
density regions around Y-pixels 40 to 90. Reassuringly a clear 
bottom-left to top-right diagonal is seen in each plot, demon- 
strating that the Fe xii and Fe xm lines are responding to the ve- 
locity shifts in unison. The effect is particularly clear for Fe xm 
/1203.83 and /1202.04. In addition to the diagonal pattern, it is 
also clear that the velocity offsets discussed above for the May 6 
data set are consistent with those for May 3, confirming the in- 
tegrity of the EIS wavelength scale over time. An interesting ef- 
fect in the velocity plots for May 3 is the loop-like structure vis- 
ible for the red-shifted pixels and most clearly seen in the Fe xm 
/I196.54-/i202.04 plot, but also partly seen in each of the other 
three plots. No explanation can be offered at this point for such 
an effect. 

In summary, the plots shown in Fig.|2T|give confidence in the 
quality of the Gaussian fits employed in the previous sections. 
In particular we see no evidence of line blending (beyond that 
already accounted for) that would lead to anomalous broadening 
or line shifts. 

9. Discussion 

Sects. |6]and|2]compared the derived densities from the four Fe xii 
and Fe xm line ratios and significant differences were seen be- 
tween the two ions. In this section we consider whether these 
differences demonstrate real physical differences between the 
Fe XII and Fe xm emitting regions, or whether atomic or instru- 
mental effects could be responsible. Firstly we consider the pairs 
of ratios from each species. 

Fexii /1196.64//1195.12 systematically yields higher densi- 
ties than /il86.88//1195.12 in both data sets. Sect. ||] identified 
a blending Fe vm line as being partly responsible, but this does 
not explain the consistent discrepancies seen in the May 3 data 
set. The /1196.64//1186.88 ratio is actually relatively insensitive 



to density (a fact used in Appendix ICjto help estimate the grating 
tilt), and Fig. |22] compares how the measured ratio values from 
the two data sets compare with the predictions from CHIANTI. 
The measured ratio values are plotted against the density de- 
rived from /1186.88//1195.12, with the predicted ratio variation 
with density from CHIANTI over-plotted. It is clearly seen that 
CHIANTI systematically under-estimates the measured ratios 
except for a small group of points in the May 3 data set. The 
shape of the gently sloping curve is reproduced quite well, ex- 
cept the high density points are a little further from the curve, 
explaining why the difference in derived densities is larger when 
the densities are high (Fig.fTTIi. The effect of the blending Fe vm 
line on /1196.64 is clearly seen through the position of the small 
group of squares in the plot which denote spatial locations where 
Fe vm is estimated to contribute more than 6 % of the intensity 
of the measured line at 196.64 A. An additional blending line 
could explain the systematic offset between the measured and 
predicted values of the /I196.64//1186.88 curve, but such a line 
would have to show very similar temperature and density be- 
haviour to /1196.64 and this is highly unlikely. 

Based on the fact that /1186.88//1195.12 gives densities closer 
to Fexm than /il96.64//1195.12 we believe that the former ra- 
tio is thus more accurate. The discrepancy between observations 
and theory revealed in Fig. |22]is then due to CHIANTI under- 
estimating the /1196.64//1195.12 ratio, however, a systematic er- 
ror in the relative calibration of EIS of around 10-15 % between 
wavelengths 186.88 A and 196.64 A can not be ruled out. 
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Fig. 22. A comparison of the measured variation of the Fexii 
/il96.64//1186.88 intensity ratio with density compared to the 
predictions from the CHIANTI database. The measured densi- 
ties are obtained from the Fexii /1186.88//1195.12 ratio. Small 
filled circles denote measurements from the May 6 data set, 
larger open circles denote measurements from the May 3 data 
set. Points from the May 6 data set where the Fe vm is predicted 
to contribute 6 % or more to the blended /1196.64 line are de- 
noted by open squares. 

The discrepancy between the densities predicted by Fexm 
A203M/A202M and A 196. 54/^202. 04 only becomes signifi- 
cant at high densities and this clearly lies in a problem with 
the /1203.82//1202.04 ratio, our reason being that the density 
variation from /1196.54//1202.04 shown in Fig. [19] tracks the 
variation from Fexii /1186.88//il95.12 very well. The densities 
from /1203.82//1202.04 increase anomalously fast compared to 
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the three other ratios when the density is above 10'*' cm"^. We 
can investigate this effect a little further by plotting the mea- 
sured /1203.82//i202. 04 ratio against the densities derived from 
/1196.54//1202.04, and comparing with the predicted ratio vari- 
ation from CHIANTI (Fig.|23ll. It is seen that above 10'° cm^^^ 
the measured ratio values are not actually far from the theoret- 
ical curve, but because the curve flattens at these densities the 
measured ratios translate to high densities using the CHIANTI 
model. A discrepancy of a similar magnitude is seen around ratio 
values of 1 .5, but this translates to only a small diff'erence in den- 
sity due to the high gradient of the theoretical curve at this loca- 
tion. A simple remedy to the high density problem is if the high- 
density limit predicted by CHIANT I is increased from t he cur- 
rent value of 4.5 to around 5.0. The lKeenan et al.l (l2007l) Fe xm 
model predicts a high-density limit of 4.7 (Fig.[T]l and so the den- 
sities derived from the May 3 data set are correspondingly lower 
in the high density region (Fig. l20l l. The high-density limit for 
/1203.82//1202.04 will be a valuable test for future atomic calcu- 
lations of Fe xm. 
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Fig. 23. A comparison of the measured variation of the Fe xni 
/1203.82//1202.04 intensity ratio with density compared to the 
predictions from the CHIANTI database. The measured densi- 
ties are obtained from the Fexiii /1196.54//1202.04 ratio. Small 
filled circles denote measurements from the May 6 data set, 
larger open circles denote measurements from the May 3 data 
set. 

Figs. [14] and [19] showed the density comparisons between 
Fe xii and Fe xm for what we consider to be the best diagnostics 
for each ion. The features in the density profiles generally match 
very well as might be expected for two ions formed at very simi- 
lar temperatures, however the absolute density values differ, with 
Fexii almost always giving higher densities than Fexm. Also 
striking in the May 6 data set (Fig.[T4b is that the range of den- 
sity variation of Fe xii is greater than that for Fe xm. 

Can the differences be due to a systematic error in the rel- 
ative calibration of the EIS instrument? By artificially varying 
the instrument sensitivity at 202 A one finds that the densities 
from Fexii ^186.88/^195. 12 and Fexm /ll96.54/^202.04 from 
the May 3 data set can be brought into line if the sensitivity at 
202 A sensitivity is increased by 35 % (i.e., the measured DN 
value at 202 A converts to a smaller line intensity, thus increas- 
ing the line ratio and the derived density). However such a large 
change would also mean that the sensitivity at the nearby wave- 



length of 203.8 A must also be significantly increased, which 
would then spoil the excellent agreement between the Fexm 
/il96.54//1202.04 and ^203.82/^202.04 ratios seen in Fig.p A 
detailed study of insensitive line ratios around 190-205 A will 
be required to confirm that the EIS relative calibration is sound, 
but we believe a discrepancy of 35 % is highly unlikely. 

Ruling out a calibration problem, we can now ask: should 
Fexii and Fexm actually yield the same density values? Even 
if the two source regions for the emitting lines do not have the 
same densities, an assumption of constant pressure would only 
lead to discrepancies of around 0. 1 dex (although the higher den- 
sity this implies for the cooler Fe xii ion is consistent with the 
measurements). A key issue to consider is whether the discrep- 
ancies could be caused through the existence of multiple den- 
sity components along the line of sight: an implicit assumption 
when generating the plots is that the plasma observed has a fixed 
density for the line of sight corresponding to each pixel. This is 
generally not the case in the corona. Even if a coronal loop, for 
example, has a fixed density at the location being considered, the 
loop is sitting in a coronal background that could have a varying 
density with height (in the case of a hydrostatic equilibrium), or 
multiple densities due to a variety of loop like structures along 
the line of sight. The measured density then actually represents 
a weighted ave rage of the different plasma components (e.g., 
iDoschekL [19841) . 

To look at this effect further, we consider the May 6 data set. 
Consider the spatial region around Y-pixel 150 where the Fexii 
and Fe xm densities become very similar. This region is seen in 
the upper panels of Figs.[T2land[T3]to be free of loops and thus 
correspond to the background corona. In the nearby pixel regions 
90-130 and 160-170 which are seen to correspond to active re- 
gion loops, the discrepancy between the two ions sharply rises 
to 0.2-0.3 dex. Is this discrepancy due to the background sub- 
traction skewing the loop density measurements? 

Fig. [24] shows the densities derived from Fexii 
^186.88//1195.12 and Fexm ^203.82/^202.04 after a back- 
ground contribution to the line intensities has been subtracted. 
For each line, the background intensity was calculated by 
averaging the intensities in Y-pixel region 132 to 143. These 
background intensities were then subtracted from the intensities 
in the 'loop' region (Y-pixels 90-120), and the densities re- 
derived. Comparing with Fig. [14] the background subtractions 
have resulted in higher densities across the loop region by about 
0.2 dex, but the discrepancy between Fe xii and Fe xm has not 
been resolved. Considering the central part of the loop region 
(Y-pixels 95-105) the average difference between Fexii and 
Fexm was 0.28 dex before background subtraction, and be- 
comes 0.27 dex after subtraction. In summary, then, background 
subtraction does not affect the Fe xii-Fe xm discrepancy in this 
case. 

The May 3 data set shows a much more complex plasma and 
so considerations of background subtraction on the density val- 
ues are more difficult. However, the density in the Y-pixel region 
60-80 is so high that the emission from this region will over- 
whelm any contribution from the background and thus we can 
speculate that it is a relatively 'pure' measurement of the region's 
density. Despite this the discrepancy between the two ions' den- 
sity values is around 0.4-0.5 dex. It is difficult to envisage any 
physical circumstances that could result in such a discrepancy. 
We also note the discrepancy over the Y-pixel region 100-175 is 
consistently at the 0.4-0.5 dex level even though there is a wide 
range of plasma structures in this region (see the upper panels of 
Figs. [TTl and [T8]). Therefore we believe these discrepancies are 
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Fig. 24. A comparison of densities derived from Fexii 
/1186.88M195.12 and Fexiii A2Q3M/A2Q2M for the loop bun- 
dle located at Y-pixels 90-120 in the May 6 data set. A back- 
ground level has been subtracted from each of the line inten- 
sities, thus leading to different densities than those shown in 
Figs.IH 

largely due to atomic physics parameters for the two ions, rather 
than physical properties of the plasma. 

A further check that can be performed on the diagnostics is 
to use the derived densities to estimate the column depth of the 
emitting region at each pixel. Comparisons of column depths 
from the two ions with the observed sizes of solar structures po- 
tentially allow a statement to be made about which ion is yield- 
ing the most accurate density values. This issue will be discussed 
in a following paper. 

10. Conclusions 

An active region was observed by Hinode/ElS on 2007 May 3 
and May 6, and coronal density values have been derived from 
emission line ratios of Fe xii and Fe xm ratios that are the most 
sensitive in the BIS wavelength bands. Methods for automati- 
cally fitting the emission lines at each pixel in the spectral im- 
ages have been described and line blending discussed. The de- 
rived densities are presented in Figs. [T2| - [T9] and the following 
conclusions are drawn. 

- In active region data sets the precision of the density mea- 
surements from the four Fe xii and Fe xiii can be as high as 
+5 % in individual spatial pixels - a considerable improve- 
ment over all previous solar spectrometers. 

- The high precision of the measurements reveals significant 
discrepancies between the four ratios, thus implying the ac- 
curacy of the measurements is only around a factor 2. 

- For Fexii, the /1196.64//1195.12 ratio always yields higher 
densities than /1186.88//il95.12 by up to 0.4 dex in log A^e- 
This is believed to be due to atomic data effects and/or in- 
strumental effects, rather than a real physical effect. 

- For Fexiii, the /1203.82//1202. 04 ratio is in very good agree- 
ment with /1196.64//1202.04 for densities less than log A^e - 
10.0. Above this density /l203.82/.i202.04 yields higher den- 
sities which is most likely due to an inaccuracy in the theo- 
retical high-density limit for the ratio. 

- Comparing Fe xn and Fe xm, the Fe xii densities are higher 
than the Fe xiii densities at almost all spatial locations, with 



the discrepancy being largest (up to 0.5 dex) at high densi- 
ties. This is most likely due to the atomic data for the ions, 
but no indication can be given here on which ion's densities 
are most accurate. 

- The May 6 data set shows a greater range of density variation 
for Fe XII than for Fe xiii which may be due to real physical 
differences for the Fe xii and Fe xm emitting regions. 

In addition to the above a number of instrumental and data 
analysis problems have been identified in this work 

- The EIS spectra are tilted relative to the CCD's axis in the 
SW band with a slope of -0.0792 A/pixel, an effect due to a 
slight misalignment of the grating relative to the CCD. This 
can have a significant effect on derived densities and must be 
corrected for if the emission lines are separated by more than 
a few angstroms. 

- It is necessary to account for the Fe xii line at 1 95 . 1 8 A when 
interpreting the intensity of the strong Fexii /1195.12 line at 
densities > 10'" cm"-'. 

- Fe VIII /1196.65 contributes to Fe xn /il96.64 in regions where 
Fe VIII emission is enhanced. The size of the contribution can 
be estimated by measuring Fe vm /II 94.66, however the cur- 
rent CHIANTI model for Fe vm appears to under-estimate 
the contribution. 

- Fe XII /1186.88 is blended with S xi /1186.84, but the contribu- 
tion of this line is found to be mostly at the level of < 10 %, 
and has only a small effect on derived Fe xii densities. 

- The rest wavelengths contained in CHIANTI for some of 
the Fexii and Fexiii lines are not correct. We find new 
wavelengths of 196.518 + 0.003 A for Fexiii /1196.54, and 
196.647 + 0.003 A for Fe xii il96.64. 

The discrepancies between the densities derived from the 
different Fe xii and Fe xm diagnostics are a serious impediment 
to using the densities for science analysis. A further study into 
the effects of the discrepancies on plasma column depths and 
filling factors will be presented in a following paper, which may 
shed light on which ion yields the most accurate density val- 
ues. We believe the discrepancies are largely due to inaccurate 
atomic parameters for one or both of the ions and so new investi- 
gations into these parameters are urgently required in order that 
scientists may make best use of the high quality EIS spectra. 
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Ion 


Wavelength/A 


Transition 




EA7cm"^ 


Fexn 


186.854 


3^23p3 2^3/2 


- 3s^3p^CP)3d ^Fs,2 


0.114 




186.887 


3,s-23p3 ^Dsi2 


-3s^3p^CP)3d^Fii2 


0.115 




195.119 


3i23p3 '*S3/2 


-3s^3p\^P)3d'*Ps,2 


0.303 




196.640 


35^3p3 ^D,i2 


-3s3p\'D)3d^Dsi2 


0.301 


Fexm 


196.540 


3s^3p^ 'D2 - 


3s^3p3d^Fi 


0.302 




202.044 




3s^3p3d^Pi 


0.084 




203.797 


3s^3p^ ^Pi - 


3s^3p3d^D2 


0.044 




203.828 


3s^3p' - 


3s^3p3d 3£)3 


0.044 



" Effective area values. 




Fig. 5. The intensity ratio of S xi /1188.64 to the blended Fe xi /1186.88 feature for two columns of spatial pixels in the May 3 (left 
panel) and May 6 (right panel) data sets. The /1188.64 has been estimated from fits to S xi /1188.68. For the May 6 data set, accurate 
fits to il 88.68 were not possible due to a dust particle on the CCD at this position. 
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Fig. 12. The upper panel shows an image from the May 6 data set obtained in the Fe xii /1195.12 line. The horizontal line indicates 
the data column chosen to derive densities. The lower panel shows densities derived from the Fexii /1186.88//1195.12 (blue) and 
/1196.64//1195.12 (red) density diagnostics. Short vertical lines are used indicating the lower and upper densities based on the Icr 
error bars on the measured line ratios. The black line shows the variation of the /1195.12 line intensity, while the dashed line shows 
the variation of the Fe vm /1186.60 line. For the latter two lines the quantity plotted is log (10^/) and log (10^^/), respectively, where 
/ is the line intensity in units erg cm"^ s ' sr ' . 




Fig. 13. The same plot as Fig. [12] except for Fe xiii. In the upper panel the image shown is from Fe xiii /1202.02. In the lower panel 
the red lines show densities derived from the /il96.54//1202.02 ratio and the blue lines shows densities from the /1203.82//1202.02 
ratio. The black line shows the variation in intensity of Fe xm /1202.04, where the plotted quantity is log (10^/). 
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Fig. 14. A comparison of densities derived from the Fexii /il86.88//il95.12 (blue lines) and Fexiii /1203.82//1202.04 (red lines) 
ratios. The black solid line shows the variation of the Fexn /i 195. 12 intensity, while the black dashed line shows the variation 
of the Fexiii /1202.04 intensity. For the latter two lines the quantity plotted is log (10^^/), where / is the line intensity in units 

erg cm"^ s"' sr"^ 
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Fig. 17. The same plot as Fig. [12] except this time for the May 3 data set. For the intensity curves from Fexn /1195.12 and Fe viii 
/II 94.66, the quantities plotted are log (10^^/) and log (10^/), respectively, where / is the line intensity in units erg cm"^ s"' sr"' . 
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Fig. 18. The same plot as Fig. [T3]except this time for the May 3 data set. For the intensity curve from Fexm /1202.04 the quantity 
plotted is log (10^^/) where / is the line intensity in units erg cm"^ s ' sr ' . 
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Fig. 19. A comparison of densities derived from the Fexii /I186.88//1195.12 (blue lines) and Fexm /1196.64//1203.82 (red lines) 
ratios in the May 3 data set. The intensity curves from Fe xii /1195.12 and Fexm /1202.04 show the quantity log (10^/) where / is the 
line intensity in units erg cm"^ s"' sr"^ 



Young et al.: High precision density measurements in the solar corona I 



19 



CD 




Fe XIII 



50 55 60 65 70 75 80 
FWHM(A,195.12)/mA 



50 55 60 65 70 75 80 
FWHM(A,195.12)/niA 




-10 10 
v(A,195.12)/km s"' 



-10 10 
v(A,195.12)/km 



CO 




65 70 75 80 85 90 
FWHM(:V195.12)/mA 



60 65 70 75 80 85 90 
FWHM(A,195.12)/mA 




10 20 
v(A,195.12)/km s 



10 10 20 
v(A,195.12)/km s ' 



-< 
E 



85 








80 


on 






■< 
E 


ID 


75 








70 












70 






CD 


65 


















65 






5 


60 








X 




60 








55 


55 








50 




50 
20 



55 60 65 70 75 80 
FWHM(A,202.04/mA 



50 55 60 65 70 75 80 
FWHM(A,202.04) / mA 




-10 10 
v(A,202.04) / km s"' 



-10 10 
v(A,202.04) / km s"' 




55 60 65 70 75 80 85 
FWHM(A,202.04) / mA 



55 60 65 70 75 80 85 
FWHM(A,202.04) / mA 




-10 
-20 
-30 
-40 
-50 



1 


1 1 


• . ■ 

:W: ■ 













10 20 
v(A,202.04) / km s" 



-10 



10 20 
v(A,202.04) / km s"' 



30 



Fig. 21. Comparisons of line widths and velocities for the density diagnostic lines of Fexii and Fexiii from the May 6 and May 3 
data sets. For Fe xii, /1195. 12 is used as the reference line, while for Fe xiii /1202.04 is used. In each plot a cross indicates the value 
of FWHM or velocity from the two lines being compared at a particular spatial pixel along the data columns discussed in Sects. |6] 
and|7] The top two rows of plots show comparisons from May 6, while the bottom two rows show comparisons from May 3. The 
two left-most columns contain the comparisons for the Fe xii lines, and the two right-most columns the comparison for the Fe xiii 
lines. 
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Appendix A: An investigation of the FeXI1 195.12 
blend 

The strong Fexii /1 195. 12 line has been known to have 
a significant blending component due to the 3s^3>ir' ^^3/2 
- 35^3p^('£ )) 3t/ ^Z) 3/2 Fexii transition since the work of 
iBinello et al.l (I2OOII) . who placed the latter transition at 
195.13 A. The wa velength was revised to 195.179 A by 
iDel Zanna & MasonI (l2b05b based on unpublished work by 
B.C. Fawcett (G. Del Zanna, private communication, 2008). In 
this section we investigate some of the properties of the fit pa- 
rameters from our two Gaussian fits. We do this by making com- 
parisons with a 1 Gaussian fit to the 195.12 A feature using the 
same wavelength pixels as marked in Fig.|6] 

First ly we compar e the lin e widths of /1 195. 12. It has been 
noted ( Doschek et al 1 I2007al) that this line is broader than 
nearby Fexii /11 93.51, and /1195.18 w as suggested as a possi- 
ble cause of this (lYoung et all lIOOTbh . In Fig. O 

we compare 

the widths of /1 193. 51 and /1 195. 12 measured from the May 3 
data set when /1195.12 is fit with a single Gaussian and when it 
is fit with a double Gaussian. The width of /1 195. 12 is clearly 
seen to be larger than that of /1193.51 with the one Gaussian fit, 
but they come into excellent agreement when a two Gaussian 
fit is used. This demonstrates that /t 195. 18 is reponsible f or the 
broader width of ^195. 12 noted bv lDoschek et al.1 (l2007al) . 

A second issue related to /J195.18 that can be investigated 
is the density sensitivity of the /I195.18//il95.12 ratio. Fig. |A.2| 
shows the predicted variation of the ratio from CHIANTI where 
/1195.18 is seen to become relatively stronger at high densities. 
Overplotted on this figure are observed ratios and densities from 
the May 6 and May 3 data sets. The densities are those obtained 
from the Fexii/il86.88//ll95.12ratio. For densities > 10'" cm"^ 
agreement is quite good, although there are a number of points 
significantly above the theoretical curve which are likely due to 
the /1195.18 intensity being over-estimated by the two Gaussian 
fits. 

At low densities, the observed ratios are all above the the- 
oretical curve although the shape of the curve is reproduced 
quite well. This may be due to a systematic over-estimate of the 
/1 195. 18 intensity in the fits due again to an unaccounted-for line. 
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Fig.A.l. A comparison of the full widths at half maxima 
(FWHM) of Gaussian fits to the Fe xii /1 1 95 . 1 2 and ^ 1 93 .5 1 lines 
in the May 3 data set. The upper panel shows the results when a 
single Gaussian fit is applied to /1 195. 12, the lower panel shows 
the results when a two Gaussian fit is applied in order to account 
for the blending /1 195. 18 Hne. 



Appendix B: The FeXII triplet 192.39, 193.51, 195.12 

The 2,s^2,p^ ^S^n - 3s^3p\^P)3d ^^^5/23/2,1/2 transitions yield 
the three Fexii transitions at 195.12, 193.51 and 192.39 A, re- 
spectively. Their theoretical intensity ratios are almost com- 
pletely insensitive to temperature and density and CHIANTI 
gives the values 1.0:0.67:0.32, res pectively. A ll three lines 
are very strong in the EIS spectra CYoung et al.L l2007 b) with 
/1195.12 the strongest and so this hne is most commonly ob- 
served with EIS. However, in active region observations /1 195. 12 
can be saturated making it difficult to use. In such circumstances 
it may be necessary to use /J192.39 or /1193.51 if they are ob- 
served with the EIS study as they are less likely to be satu- 
rated. In this section the relative intensities as measured from 
the May 3 and May 6 data sets are presented and compared with 
the CHIANTI predictions. 

The ratios of the three lines will clearly be affected by blend- 
ing with other species, and the contribution of /1 1 95 . 1 8 to the fea- 
ture at 195.12 A was highlighted in Appendix lAl which demon- 
strated that /1195.18 is responsible for the enhanced broadening 
of the 195.12 A feature when it is fit with a single Gaussian. 



Inspection of the line widths of /1192.39 and /1193.51 reveals no 
anomalies consistent with a lack of blending for these lines. 

Fig. 10 shows the il92.39//il95.12 and /1193.51M195.12 
intensity ratios from the May 3 and May 6 data sets (lower and 
upper panels, respectively). Despite the large variation in density 
and intensity for the data sets, the ratios show excellent agree- 
ment with the CHIANTI theoretical values. The only discrep- 
ancy is that the /1193.51//il95.12 ratio in the May 6 data set is 
noticeably above the CHIANTI prediction and the May 3 val- 
ues, and no explanation can be given for this effect. Another 
effect that can be seen in both ratios is that the measurements 
rise slightly as the Y-pixel increases. This can be explained if 
the instrument sensitivity is slightly different towards the top of 
the EIS slit than at the bottom, becoming more sensitive at short 
wav elengths at the top. The laboratory measurements performed 
by Seelv et aTl ( 1200 4) showed the shape of the multilayer coat- 
ing reflectivity varying with the position of the light source in the 
aperture, and thus the Fe xii ratios could be revealing this effect. 

In summary, though, the /1192.39//1195.12 and 
/1193.51//1195.12 ratios show excellent agreement with 
theoretical predictions and thus either /1192.39 or /1193.51 can 
be used in place of /1195. 12 when forming density diagnostics. 
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Fig. A.2. A comparison of the measured variation of the Fe xii 
/1195.18//1195.12 intensity ratio with density compared to the 
predictions from the CHIANTI database. The measured densi- 
ties are obtained from the Fexii /1186.88//1195.12 ratio. Small 
filled circles denote measurements from the May 6 data set, 
while larger, open circles denote measurements from the May 3 
data set. 



Appendix C: Correcting for spatial offsets between 
emission lines 

The spectral image formed on the SW CCD is marginally tilted 
relative to the CCD axes such that the same spatial feature is 
observed to be slightly lower on the CCD in long wavelength 
lines compared to short wavelength lines. The effect is due to a 
misaligment of the grating relative to the CCD and we refer to 
it here as the "grating tilt". The best way to measure the tilt is 
to identify emission lines formed at the same temperature and 
which are density insensitive relative to each other. If a com- 
pact spatial feature is observed in such lines, then the offset can 
be estimated by co-aligning the images in the two lines. For the 
EIS SW band, however, there are few emission line pairs that 
are insensitive and span a wide enough wavelength range to be 
useful. In addition, most of the lines are formed in the corona 
where spatial features are more diff'use than in the transition re- 
gion, making alignment more difficult. A further complication 
is that the tilt is small (a one spatial pixel offset over 12.6 A is 
found here), meaning for most line pair comparisons the offset 
is on the sub-pixel level. 

For the present work the focus is on the density sensitive 
lines of Fexii and Fexm which are found between 186.8 and 
203.8 A. We thus seek to identify the grating tilt over this re- 
gion. We make use of the May 3 and May 6 data sets that have 
been used for the density analysis, and four line pairs are con- 
sidered. Fe VIII images oft en show compact fea tures that are very 
suitable for co-alignment (I Young et al.Ll200'7ah . and the /1185.21 
and /1194.66 lines are used here. They have a wide wavelength 
separation and only a small density sensitivity, however /1185.21 
is blended with a Ni xvi line and so care has to be taken to avoid 
comparing spatial features where there is clearly hot emission. 

The remaining three ratios are all from Fexii. Fexii 
/1186.88//1196. 64 involves two of the lines from the density anal- 
ysis. CHIANTI predicts the intensity ratio to vary from 0.21 
to 0.27 over the density range 10**-10'^ cm"^ and so the ra- 
tio is relatively insensitive to density. Fexii /1192.39//1195.12 
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Fig. B.l. Comparison of measured Fexii /1193.51//1195.12 and 
/1192.39//il95.12 intensity ratios against the predictions from 
CHIANTI from the May 6 (upper panel) and May 3 (lower 
panel) data sets. The solid horizontal lines in each panel show 
the line ratio values predicted from CHIANTI, and crosses de- 
note the measured ratio values. The average of the measured val- 
ues is shown for each of the ratios and data sets. 

involves two very strong lines that are insensitive to density 
(Appendix IBTi. however the wavelength separation is small. The 
last ratio is /i203.73//ll95.12 which, although it shows signifi- 
cant density sensitivity, was deemed the best means of estimat- 
ing the grating tilt out to the spectral feature at 203.8 A. /1203.73 
is measured through a multi-Gaussian fit to Fexm /1203.82 
(Sect.l53Tl. 

The method for finding the spatial offset between two lines is 
as follows. Each of the lines was fit with a Gaussian as described 
in Sect.|5j yielding intensity images in each of the lines. A com- 
pact feature was then identified in particular image columns, and 
the intensity profiles along these columns extracted. One of the 
two columns was treated as a reference and, for the other, the 
intensity values were adjusted to simulate moving the image up 
and down on the detector. The intensity for this comparison col- 
umn was re-calculated as 



Aiew(0 



(\ - a)l{i) + al(i + Y) a > 
(1-|q'|)/(0 + |Q'|/(/-1) Q'<0 



(C.l) 



where a is varied from -1 to +1 in 0.1 intervals. There are thus 
21 columns to be compared with the reference column. From 
the reference column, the portion containing the compact spa- 
tial feature was extracted and compared with the 21 compari- 
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son columns to determine which gave the best spatial match. 
This was done by first normalising the extracted columns to the 
maximum value in the extracted region, and then calculating the 
quantity: 



(C.2) 



where 7*^^ is the normalised intensity of the reference wave- 
length, and is the normalised comparison intensity 
(Eq. IC.lb . The best match to the reference column was deter- 
mined to be the comparison column with the lowest value. 

The May 6 data set had limited value for estimating the spa- 
tial offsets as Fe xn does not demonstrate any compact spatial 
features (e.g., the upper panel of Fig. [T7] i. However, there is a 
bright, compact feature seen in Fe vm (lower panel of Fig. [TTl l. 
The offset from the /II 85.21, /II 94.66 comparison was found to 
be 0.7 pixels. 

The May 3 data set showed a number of compact features at 
different temperatures however the Fe viii comparison was com- 
promised by the blending of /II 85.21 and so the offset from the 
May 6 dataset was preferred. Comparing several features, the 
offset between Fexii /1186.88 and /1196.64 was found to vary 
from 0.7 to 0.8 pixels - we assume 0.75 here. For /1195.12 and 
/1203.72 the offset was 0.8 pixels, and for /1192.39 and A\95.\2 
it was 0.2 pixels. 

Assuming that the grating tilt is linear, the offset for a given 
wavelength is written as 



y = m(/l- 195.12) 



(C.3) 



where the offset is assumed to be zero at 195.12 A. Computing 
the tilt from the four ratios and taking the average, we find m - 
-0.0792 pixels/A. 

To illustrate the consequence of the grating tilt on the derived 
densities, we consider a small brightening that appears at X-pixel 
25 and Y-pixel 98 in the May 3 data set (see Fig. Fig. ICT] 
compares the densities derived from the two Fe xii density diag- 
nostics in the vicinity of this brightening. The top panel shows 
the densities when no offsets are applied, the bottom panel shows 
the densities when the offsets have been applied. It is clearly seen 
that the peaks in density around pixels 96-100 do not coincide 
in the top panel, but do coincide in the bottom panel. The dis- 
crepancy in the top panel is striking and greater than would be 
expected simply from the sub-pixel spatial offsets derived for 
the lines. This is because taking the ratios of lines amplifies the 
effect of the intensity offsets. 

Finally we note that there is some uncertainty in deriving the 
grating tilt due to the difficulty of co-aligning data at the sub- 
pixel level. This thus represents a possible systematic error in 
the derived density measurements that will be greatest for those 
lines with a large wavelength separation. It may affect density 
comparisons between ratios for small spatial features, but not 
the broad trends highlighted in Sect.|9] 
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Fig. C.l. A comparison of densities obtained from the Fexii 
/i 1 86 . 8 8//1 1 95 . 1 2 and /1 1 96 .64//i 195.12 ratios measured from the 
May 3 data set. A section of a data column from X-pixel 25 (see 
upper panel of Fig. [T7]for reference) is shown where is found a 
compact high density feature. The upper panel shows the densi- 
ties obtained when the grating tilt is not corrected for, and the 
lower panel the results when the tilt is corrected for Note that 
the discrepancy between the density peaks in the upper panel is 
resolved in the lower panel. 
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Fig. 15. Exactly analogous plot to Fig.[T3] except the densities are derived using the lKeenan et alJ (l2007h Fe xiii atomic model. 
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Fig. 20. Exactly analogous plot to Fig. [18] except the Fe xiii densities are derived using the lKeenan et alj (|2007|) atomic model. 



